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Background: Potential energy surfaces (PES’s) of actinide nuclei are characterized by a two-humped barrier 
structure. At large deformations beyond the second barrier the occurrence of a third one was predicted by 
macroscopic-microscopic model calculations in the 1970s, but contradictory results were later reported by number 
of studies that used different methods. 


Purpose: Triple-humped barriers in actinide nuclei are investigated in the framework of covariant density func- 
tional theory (CDFT). 


Methods: Calculations are performed using the multidimensionally-constrained relativistic mean field (MDC- 
RMF) model, with the nonlinear point-coupling functional PC-PK1 and the density-dependent meson exchange 
functional DD-ME2 in the particle-hole channel. Pairing correlations are treated in the BCS approximation with 
a separable pairing force of finite range. 


Results: Two-dimensional PES’s of 726:278:299232Th and 732:234,236.238 U are mapped and the third minima on 
these surfaces are located. Then one-dimensional potential energy curves along the fission path are analyzed in 
detail and the energies of the second barrier, the third minimum, and the third barrier are determined. The 
functional DD-ME2 predicts the occurrence of a third barrier in all Th nuclei and ?°°U. The third minima in 
230,232Th are very shallow, whereas those in 72°??8Th and 7°°U are quite prominent. With the functional PC- 
PK1 a third barrier is found only in 7?°??5:?2°Th. Single-nucleon levels around the Fermi surface are analyzed in 
226Th, and it is found that the formation of the third minimum is mainly due to the Z = 90 proton energy gap 
at B20 x 1.5 and B30 = 0.7. 


Conclusions: The possible occurrence of a third barrier on the PES’s of actinide nuclei depends on the effective 
interaction used in multidimensional CDFT calculations. More pronounced minima are predicted by the DD-ME2 
functional, as compared to the functional PC-PK1. The depth of the third well in Th isotopes decreases with 
increasing neutron number. The origin of the third minimum is due to the proton Z = 90 shell gap at relevant 
deformations. 


PACS numbers: 21.60.Jz, 24.75.+i, 25.85.-w, 27.90.+b 


I. INTRODUCTION measurements for ?30-233Th and 23"U supported the ex- 


istence of shallow third minima on the PES’s of these nu- 


The potential energy surface (PES) of a fissile nucleus le! [9-15]. The deduced values for the depth of the third 


plays a crucial role in fission studies. Observables such 
as the fragment mass distribution and the fission half-life 
are closely related to properties of the PES of a com- 
pound nucleus. The PES’s of actinide nuclei are char- 
acterized by a two-humped barrier structure along the 
static fission path, which has extensively been studied 
both experimentally and using a variety of theoretical 
models. We refer the reader, for instance, to Refs. [1—4] 
and references therein. Macroscopic-microscopic model 
calculations predicted, already in the 1970s [5-7], the 
occurrence of shallow third minima beyond the second 
barrier. These minima were used to explain the tho- 
rium anomaly [6-8]. High resolution fission cross section 
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well are only a few hundred keV. In Ref. [16] a model 
was developed to describe fission in light actinides, and 
to consider transmission through a triple-humped fission 
barrier with absorption. Using this model the complex 
resonance structure in the first-chance neutron-induced 
fission cross sections of ?32Th and 2%!Pa was reproduced, 
and shallow third minima with a depth of less than 1 
MeV were obtained. 


The PES’s of nuclei in the Ra-Th region were com- 
puted using the macroscopic-microscopic model with a 
modified oscillator potential, and in many nuclei a third 
minimum was found at very large quadrupole defor- 
mation [17]. It was concluded that the depth of the 
third minimum could be at most 1.5 MeV. This model, 
with shell corrections calculated by adopting an axially- 
deformed Woods-Saxon potential, was later used to sys- 
tematically study the PES’s of even-even Rn, Ra, Th, 
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and U isotopes [18]. Very deep minima, or even two hy- 
perdeformed minima, were predicted in many of these 
nuclei. In some cases the depth of the third minimum 
could be as large as 4 MeV. 


A series of experiments were performed to find evi- 
dence for hyperdeformed states in U isotopes [19-23]. 
The deduced values for the excitation energies of the 
third minima and the third barriers were Erm = 3-4 MeV 
and Brim ~ 6 MeV, respectively. Thus, the depth of the 
third well could be 2-3 MeV. Such deep third minima in 
U isotopes agree with the predictions of Ref. [18], but dif- 
fer from the experimental results of 7°°-?33Th and 287U 
[9-15], and from the theoretical predictions reported in 
Refs. [5-7, 17]. 


To verify the predictions for third minima in ac- 
tinides in the macroscopic-microscopic model based on 
the Woods-Saxon potential, calculations with additional 
shape degrees of freedom were preformed in Refs. [24, 25]. 
It was found that the third barrier could be lowered 
substantially by including the (6, deformation. Conse- 
quently, third minima disappeared in many nuclei, ex- 
cept for ?30:-232Th, in which only a shallow third minimum 
with a depth of less than 400 keV was found. Further- 
more, in Ref. [26] an analysis of PES’s computed with 
the finite-range liquid-drop model [27] revealed that only 
few nuclei accessible to experiment exhibit third minima 
in their PES’s and the depth of the third well is less than 
1 MeV for the light Th and U isotopes. 


In addition to the macroscopic-microscopic model, a 
number of self-consistent approaches have also been used 
to investigate PES’s of deformed nuclei. Deformation- 
constrained Hartree-Fock or Hartree-Fock-Bogoliubov 
calculations with Skyrme forces [28-30] and the Gogny 
force [31-33] did not exhibit deep third minima in ac- 
tinide nuclei. In Ref. [34] a shallow third minimum with a 
depth of 1-2 MeV in ?3?Th was found in axially-deformed 
relativistic mean-field calculations with the effective in- 
teractions PL-40, NL1, and NL-SH. 


In the present study we examine the occurrence and 
properties of third minima on the PES’s of light ac- 
tinides, using the multidimensionally-constrained rela- 
tivistic mean-field (MDC-RMF) model. Third minima, 
if they exist, are located at very large quadrupole and 
octupole deformations. It is probable that some intruder 
high-N single-particle states play an important role in 
the formation of these minima. Therefore, in addition to 
the PES’s, we analyze in detail the single-particle level 
structure at hyperdeformation, and study the origin of 
possible third minima. 


The paper is organized as follows. The MDC-RMF 
model is introduced in Sec. II. In Sec. IIT we present nu- 
merical details and the results for the PES’s and third 
barriers in light actinide nuclei. A summary and conclu- 
sions are given in Sec. IV. 


Il. THE MDC-RMF APPROACH 


Relativistic mean-field (RMF) models have been ap- 
plied very successfully in studies of a variety of nuclear 
structure phenomena. In the standard representation of 
RMF models a nucleus is described as a system of Dirac 
nucleons coupled to exchange mesons through an effective 
Lagrangian [35-43]. At the energy scale characteristic for 
nuclear binding and low-lying excited states, the meson 
exchange is just a convenient representation of the effec- 
tive nuclear interaction, and can be replaced by the corre- 
sponding local four-point (contact) interactions between 
nucleons [44, 45]. A quantitative treatment of nuclear 
matter and finite nuclei necessitates a medium depen- 
dence of effective mean-field interactions that takes into 
account higher-order many-body effects. The medium 
dependence can be introduced either by including non- 
linear terms in the Lagrangian [46-48] or by assuming an 
explicit density dependence for the meson-nucleon cou- 
plings [49, 50]. Therefore, the framework of covariant 
density functionals can be formulated in one of the fol- 
lowing four representations for the effective nuclear in- 
teraction: meson exchange (ME) or point-coupling (PC) 
nucleon interactions combined with nonlinear (NL) or 
density-dependent (DD) vertex functions. 

In a self-consistent mean-field approach that allows 
breaking both axial quadrupole and reflection symme- 
tries, the MDC-RMF model has recently been developed 
and applied to studies of PES’s and fission barriers of 
actinides [1, 51-54], shapes of hypernuclei [55, 56], and 
nonaxial-octupole Y32 correlations in N = 150 isotones 
[57]. The nuclear shape is parameterized by the defor- 
mation parameters 


At 


Pru = SIR 


(Qan); (1) 


where Qanu = eee is the mass multipole operator. The 
shape is invariant under the exchange of the x and y 
axes and, thus, all deformations §),, with even u can be 
included simultaneously. The deformed RMF equations 
are solved by an expansion in the axially-deformed har- 
monic oscillator (ADHO) basis [58, 59]. For details of 
the MDC-RMF model, we refer the reader to Ref. [1]. 


II. RESULTS AND DISCUSSION 
A. Implementation and numerical details 


In the present study we employ two relativistic den- 
sity functionals: PC-PK1 [60] with point-coupling effec- 
tive interactions that include nonlinear terms in the self- 
energies and the meson-exchange functional DD-ME2 
[61] with density-dependent meson-nucleon couplings. 
Pairing correlations are taken into account in the BCS ap- 
proximation with a separable pairing force of finite range 
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FIG. 1. (Color online) Energy curve of ??°Th computed with 
the MDC-RMF model using the functional PC-PK1. When 
axial symmetry is imposed, the energy denoted by the dash- 
dotted (black) curve is obtained [axial symmetric (AS) and re- 
flection asymmetric (RA)], whereas triaxial calculations yield 
the solid (blue) curve [non-axial (NA) and reflection asym- 
metric (RA)]. The ADHO basis contains Ny = 16 oscillator 
shells. 


[62-64] 


JŘ- RPP HSE, (2) 


where g is the pairing strength and R and ř are the 
center-of-mass and relative coordinates of the two nucle- 
ons, respectively. P(r) denotes the Gaussian function: 


1 —r? jaa? 

P(r) = (4na? 372° a (3) 
where a is the effective range of the pairing force. The 
two parameters g = go = 728 MeV-fm? and a = 0.644 
fm [62, 63] have been adjusted to reproduce the density 
dependence of the pairing gap at the Fermi surface in 
symmetric nuclear matter and calculated with the Gogny 
force D1S. 

Nonaxial shapes are crucial for determining the height 
of both the inner [65, 66] and outer barriers [1, 51-54]. 
What role do they play at the third barrier? In Fig. 1 
we display the energy curve of ?2°Th computed with and 
without the inclusion of triaxial deformations. This cal- 
culation is performed in an ADHO basis truncated to 
Np = 16 oscillator shells (see Ref. [1] for the explanation 
of the truncation). As shown in the figure, triaxial defor- 
mations lower the second barrier considerably. Beyond 
the second saddle point the influence of nonaxial defor- 
mations on the binding energy gradually vanishes as (29 
increases, and these shapes appear not to be important 
at the third minimum and the third barrier. 

In Ref. [1] the convergence of the total energy with 
respect to the size of the ADHO basis has been verified 
along the fission path up to the second fission barrier. 
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FIG. 2. (Color online) The axially symmetric energy curve of 
226Th computed with the MDC-RMF model, using the func- 
tional PC-PK1. Results obtained with different truncations of 
the ADHO basis, i.e., with Ny = 16, 18, 20, and 22 shells are 
plotted by the dash-dotted, dotted, dashed, and solid curves, 
respectively. 


Since here we consider a region with even larger defor- 
mation, we will extend this check and examine the basis 
truncation up to the third fission barrier. The energy 
curve of ??6Th is calculated up to B29 = 2.7 assuming ax- 
ial symmetry and employing ADHO bases with Ny = 16, 
18, 20, and 22 shells. In Ref. [1] it has been shown that 
a basis with Nf = 20 produces very accurate results for 
the height of the second barrier, around which the triaxial 
deformation is also important. In Fig. 1 one notices that 
nonaxial shapes do not influence the height of the third 
barrier, and so the axial symmetry is imposed and we 
mainly focus on the deformation region beyond the sec- 
ond barrier. Figure 2 shows that the height of the third 
barrier is lowered when Ny increases from 16 to 20. The 
results obtained with Ns = 20 and 22 are almost identi- 
cal, and this means that Ns = 20 should present an suffi- 
cient choice. We notice that for the hyperdeformed min- 
imum around {29 ~ 1.6, results obtained with Ny = 18, 
20, and 22 are difficult to differentiate. 


B. Two-dimensional PES’s 


We consider light even-even actinides 732:294236.238 


and 276,228,230,232Th_ Figure 3 displays the self-consistent 
MDC-RMF energy surfaces in the (820, 830) plane, calcu- 
lated with the relativistic functionals PC-PK1 and DD- 
ME2. The deformation parameters are restricted to the 
range 329 ~ (1.0,3.0) and $39 ~ (0.3,1.5), in which the 
second barrier together with the third minimum and the 
third barrier, if they exist, are located. Note that the 
contour interval for the PES’s plotted in Fig. 3 is 0.5 
MeV. 

In Fig. 3(a) one notices that for 
PES’s do not display a third minimum when calculated 


232,234,236,238 77 the 
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FIG. 3. (Color online) Self-consistent MDC-RMF energy surfaces of U and Th isotopes in the (620, 830) plane, calculated 
with the relativistic functionals PC-PK1 (a) and DD-ME2 (b). For each nucleus the energy is normalized with respect to 
the binding energy of the absolute minimum. The contours join points on the surface with the same energy and the energy 
difference between neighboring contours is 0.5 MeV. The calculation has been performed in the ADHO basis with Ny = 16 
shells. 
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with PC-PK1. For these four U isotopes the second sad- 
dle point is located at B29 ~ 1.2-1.3 and B39 ~ 0.3- 
0.4. With the increase of Bao the total binding energy 
decreases monotonically along the lowest fission path. 
Different results are obtained with DD-ME2, as shown 
in Fig. 3(b). With the exception of SU, we find a 
third minimum on the PES’s of ?32:234.238U, The min- 
imum is rather shallow for ?32?U and ?34U and very pro- 
nounced for ?°°U. These results are very similar to those 
obtained in recent calculations based on the macroscopic- 
microscopic model [24—26, 67] and the Skyrme Hartree- 
Fock-Bogoliubov model [30]. 

In the case of Th isotopes, both PC-PK1 and DD-ME2 
predict. a pronounced third minimum for ??°Th with a 
depth of 2-3 MeV. As the neutron number N increases, 
both the energy of the third minimum and the height 
of the third barrier decrease and the depth of the third 
well is also reduced. This trend has also been predicted 
in the macroscopic-microscopic model calculations [26]. 
For 72°Th only a shallow minimum of depth less than 1 
MeV occurs. The third minimum completely disappears 
for ?°?Th when calculated with PC-PK1 and it is very 
shallow with the functional DD-ME2. 

Figure 3 clearly shows that the functional DD-ME2 
(lower panel) predicts more pronounced third minima 
when compared to those obtained with PC-PK1 (upper 
panel). In addition to the lowest fission path on which we 
focus in this work, there are other possible paths along 
which there are more shallow minima and saddle points. 
It would be worthwhile to investigate these fine struc- 
tures in a future study. 


C. Pairing strength 


In general the height of fission barriers is rather sen- 
sitive to the strength of the pairing interaction [68]. As 
explained above, the parameters of the separable pairing 
force of finite range that we use in this work were origi- 
nally adjusted to reproduce the pairing gap at the Fermi 
surface in symmetric nuclear matter and calculated with 
the Gogny force D1S. A number of studies based on the 
relativistic Hartree-Bogoliubov model have shown that 
this pairing force can be used to calculate structure prop- 
erties with a success; but in some other cases, the pairing 
strength needs to be fine-tuned [69, 70]. 

In the present work, since pairing correlations are 
treated in the BCS approximation, it is necessary to ver- 
ify whether the strength of the pairing force is adequate 
for the mass region under consideration. We have there- 
fore calculated the odd-even mass differences for the Th 
isotopes: 


Ay (Z,N) = = [E(Z, N +1)+ E(Z, N —1)-2E(Z,N)], 


Nl rR N| =e 


A,(Z,N) = =(E(Z+1,N)+ E(Z —- 1, N) -2E(Z,N)], 


using different pairing strengths: g/go = 1.0, 1.1 and 


—@— g/g, =1.0 —A— g/g =1.1 —*— g/g =1.2 Exp. 


FIG. 4. (Color online) The odd-even differences of binding en- 
ergies of Th isotopes computed with the MDC-RMF model, 
using the functional PC-PK1. Results obtained using differ- 
ent pairing strength parameters: g/go = 1.0, 1.1, and 1.2, 
where go = 728 MeV-fm°, are plotted as dots, triangles, and 
stars respectively. The squares denote the experimental val- 
ues obtained from Ref. [71]. The calculation has been per- 
formed in the ADHO basis with Ny = 20 shells. 


1.2, where go = 728 MeV-fm?. The results are shown 
in Figs. 4 and 5. For the case of PC-PK1, when com- 
pared to experimental values obtained from Ref. [71], 
one notices that the calculation with g/go = 1.0 under- 
estimates both proton and neutron odd-even mass dif- 
ferences considerably. While the results obtained with 
g/go = 1.1 nicely reproduce the empirical proton odd- 
even mass differences, the experimental neutron gaps are 
located between the values obtained using g/go = 1.1 and 
g/go = 1.2. For the case of DD-ME2, the results obtained 
with both g/go = 1.0 and g/go = 1.1 underestimate pro- 
ton and neutron odd-even mass differences. The results 
obtained with g/go = 1.2 reproduce the empirical proton 
odd-even mass differences quite well, and slightly under- 
estimate the neutron odd-even mass differences. Since 
we do not wish to introduce additional model parame- 
ters by considering different pairing strengths for protons 
and neutrons, in the remainder of this study we will use 
g/go = 1.1 both for protons and neutrons for PC-PK1 
and g/go = 1.2 for DD-ME2. 


Multidimensional self-consistent 
constrained calculations are very 


deformation- 
time-consuming 
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FIG. 5. (Color online) The odd-even differences of binding en- 
ergies of Th isotopes computed with the MDC-RMF model, 
using the functional DD-ME2. Results obtained using differ- 
ent pairing strength parameters: g/go = 1.0, 1.1, and 1.2, 
where go = 728 MeV-fm?, are plotted as dots, triangles, and 
stars respectively. The squares denote the experimental val- 
ues obtained from Ref. [71]. The calculation has been per- 
formed in the ADHO basis with Ny = 20 shells. 


but, on the other hand, to accurately locate the saddle 
points it is necessary to consider several shape degrees 
of freedom. Therefore, to locate the fission paths, 
in the present study we adopt the approach used in 
Refs. [1, 51]: (1) From the two-dimensional energy 
surfaces on the (29, 839) plane shown in Fig. 3, and 
calculated with Ny = 16 and g/go = 1.0, one approxi- 
mately identifies the lowest fission path; (2) In a second 
step, one-dimensional constrained calculations with 
Ny = 20 and g/go = 1.1 for PC-PK1 and g/go = 1.2 
for DD-ME2 are performed along this approximate 
fission path for nuclei with apparent third minima. In 
these calculations additional nonaxial shapes are allowed 
around the second barrier. 


D. The third minima and barriers 


By examining the two-dimensional PES’s shown in 
Fig. 3, one notices that third minima and barriers appear 
for 226,228,230,232TH and ?32:234.238U, For these nuclei the 
energy curves along the lowest static fission path, calcu- 


TABLE I. Excitation energies (in MeV) of the second sad- 
dle point Bu, the third minimum Erm, and the third sad- 
dle point Bur, with respect to the deformed mean-field 
equilibrium state for 276225230,282-Th and ?32234238U. ob- 
tained from MDC-RMF calculations. AE = Bry — Erm de- 
notes the depth of the third well relative to the third bar- 
rier. The empirical values (denoted by “Emp”) are from 
Refs. [11, 15, 20, 22, 23, 73]. 


Nucleus Parameters Bu Ey Bir AE 
26TH DD-ME2 8.76 7.37 9.31 1.94 
PC-PK1 7.94 544 6.73 1.29 
228 Th DD-ME2 8.16 669 7.82 1.13 
PC-PK1 7.19 4.72 5.50 0.78 
230TH DD-ME2 7.84 5.97 6.60 0.63 
PC-PK1 6.56 4.01 445 0.44 
Emp [73 6.80 
Emp [11 5.75 5.55 6.45 0.90 
232° Th DD-ME2 7.53 542 5.92 0.50 
Emp [73 6.70 
Emp [15 6.50 5.40 5.70 0.30 
23277 DD-ME2 7.25 
Emp [73 5.40 
Emp [22 4.91 3.20 6.02 2.82 
2347] DD-ME2 7.01 
Emp [73 5.50 
Emp [20 31 
2387 DD-ME2 7.70 3.70 481 1.11 
Emp [73 5.50 
Emp [23 5.6 3.6 5.6 2.0 


lated in the ADHO basis with Ny = 20 shells and with 
the enhanced pairing strength g/go = 1.1 for PC-PK1 
and g/go = 1.2 for DD-ME2, are shown in Fig. 6. In 
the vicinity of the second saddle point, MDC-RMF cal- 
culations are performed with and without the inclusion of 
triaxial quadrupole shapes and both results are displayed 
for comparison. As already reported in Ref. [51], the in- 
clusion of triaxial configurations in addition to the ax- 
ial octupole deformation, modifies the shape and height 
of the second fission barrier. We should note that the 
effect of triaxiality on the second barrier sensitively de- 
pends on the strength of the pairing interaction. For in- 
stance, for the case g/go = 1.0 (cf. Fig. 1), the inclusion 
of triaxial configurations lowers the second fission bar- 
rier of ?2©Th by 0.72 MeV, whereas for g/go = 1.1 (cf. 
Fig. 6) the second barrier of ?7°Th is only lowered by 
0.33 MeV. By further increasing the pairing strength the 
influence of nonaxial deformations on the binding energy 
near the second fission barrier may be further reduced 
[72]. The energies of the third minima, the heights of the 
second and third barriers, and the depths of the third 
wells AF = By — Erm, together with empirical parame- 
ters [11, 15, 20, 22, 23, 73] are included in Table I. 

For ?2°Th the functional DD-ME2 (PC-PK1) predicts 
the third minimum at 7.37 (5.44) MeV above the de- 
formed ground state. With DD-ME2 the third barrier is 
slightly higher than the second one, whereas the opposite 
is obtained with PC-PK1. The depth of the third well 


PC-PK1 | 


axial sym. 
triaxial 


FIG. 6. (Color online) MDC-RMF energy curves of 276:2?8:739282-Th and ?32234238U, The calculation has been performed in 
the ADHO basis with N+ = 20 shells and the pairing strength parameter g/go = 1.1 for PC-PK1 and g/go = 1.2 for DD-ME2. 
For each nucleus the energy is normalized with respect to the binding energy at the absolute minimum. 


computed with DD-ME2 and PC-PK1 is 1.94 MeV and 
1.29 MeV, respectively. In the case of ??STh the third 
minimum is shallower: 1.13 MeV for DD-ME2 and 0.78 
MeV for PC-PK1. Among the nuclei considered here, 
these two isotopes exhibit the most pronounced third 
minima. The fission barrier parameters for ?3°Th de- 
duced in Ref. [11] are Erm = 5.55 MeV and Bim = 6.45 
MeV. MDC-RMF calculations with the functional DD- 
ME2 reproduce these values, whereas they are underesti- 
mated by about 2 MeV by PC-PK1. For ?°?Th the third 
minimum appears on the PES calculated with DD-ME2 
but not on the one obtained with PC-PK1. The electron- 
induced fission cross section measurement for 7°?Th indi- 
cates a shallow third minimum with a depth of about 0.30 
MeV [15]. The theoretical result obtained with DD-ME2 
(0.50 MeV) is consistent with this empirical value. In 
general, Table I shows that the depth of the third min- 
imum in the Th isotopes decreases with increasing the 


neutron number. 

For the U isotopes MDC-RMF calculations with the 
functional DD-ME2 predict the existence of the third 
minimum in ?38U. The depth of the third well is 1.11 
MeV, which is smaller than the empirical value 2.0 
MeV [23]. The reason for that is shown in Table I, 
where one notices that the calculated energy of the third 
minimum (3.70 MeV) is close to the empirical value (3.6 
MeV), but the theoretical height of the third barrier (4.81 
MeV) is much smaller than the empirical value 5.6 MeV. 


E. Single-nucleon level structure 


The appearance of a hyperdeformed minimum must 
have its origin in the single particle level structure. 
Since ??6Th displays the most pronounced third mini- 
mum among the nuclei considered in this work, we will 
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FIG. 7. (Color online) The neutron (upper panel) and proton (lower panel) single-particle levels of ??°Th near the Fermi surface 
along the static fission path. For 820 < 0.6, only reflection-symmetric deformations are considered and the red (blue) curves 
represent positive (negative) parity states. When 629 > 0.6, the octupole deformation (30 has a non-vanishing value and parity 
is not a good quantum number. The dash-dotted (green) curves denote the Fermi energy, and the red (blue) symbols in the 
lower panel are used to guide the eye. The MDC-RMF calculation has been performed with the functional DD-ME2. 


analyze the neutron and proton deformed single-particle 
levels of this isotope. 


In Fig. 7 we display the neutron and proton de- 
formed single-particle levels of ??°Th near the Fermi sur- 
face along the static fission path, as functions of the 
quadrupole deformation (629. The levels are obtained in 
a MDC-RMF calculation using the functional DD-ME2. 
The quadrupole deformation of the superdeformed min- 
imum of ?76Th is predicted at B29 ~ 0.6. For Ba9 < 0.6, 
octupole deformations are not considered, and thus par- 
ity is conserved. When {29 > 0.6, the octupole deforma- 
tion {39 is nonzero and the parity cannot be considered 
as a good quantum number for large deformations. Fur- 
thermore, around the second saddle point triaxial defor- 
mations also play a role and the third component of the 
angular momentum is not conserved either. This results 
in a very complex single-particle level scheme around the 
second barrier. Since in this study we are mainly inter- 


ested in single-particle levels in the region of the third 
minimum and the third barrier, in Fig. 7 we only plot 
results obtained by imposing axial symmetry. 


The hyperdeformed minimum of ?2°Th is located at 
Boo ~ 1.5 and 639 ~ 0.7. By examining the neutron 
single-particle levels around $29 = 1.5, in the upper panel 
of Fig. 7 one notices a region of low level density near the 
Fermi surface, even though the energy gap is not large. 
For protons, as shown in the lower panel of Fig. 7, a 
large energy gap is clearly visible at Z = 90 in the region 
B20 © 1.5. Therefore, the formation of the third mini- 
mum on the PES of ??°Th is mainly caused by the large 
proton gap at the Fermi surface. Many single-particle 
states around the proton Fermi level are involved in the 
formation of the energy gap at Z = 90. These states 
are dotted with red (blue) symbols in the lower panel of 
Fig. 7, and labeled with Q, i.e., the third component of 
the angular momentum. 
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IV. SUMMARY 


We have analyzed the energy surfaces of light even- 
even Th and U isotopes using the multidimensionally- 
constrained relativistic mean field (MDC-RMF) ap- 
proach. Calculations have been performed with two rela- 
tivistic density functionals: PC-PK1 with point-coupling 
effective interactions that include higher order terms, and 
the meson-exchange functional DD-ME2 with density- 
dependent meson-nucleon couplings. Pairing correlations 
are taken into account in the BCS approximation with a 
separable pairing force of finite range. 


In a first step we have examined the two-dimensional 
PES’s of ?26-228,230,232Th and 232,234,236,238 7] and located 
the third minima on the energy maps. By analyzing 
the resulting energy curves along the lowest static fis- 
sion path, the energies of the second barrier, the third 
minimum, and the third barrier have been determined. 
In calculations with the functional DD-ME2, a third po- 
tential barrier is predicted in all Th nuclei and ?38U. The 
third well in ?3%232Th is very shallow with a depth of less 
than 1 MeV, whereas the third well in ?26-228Th and ?°°U 
is rather deep. The functional PC-PK1 predicts a third 
barrier only in 726:228.230Th, Therefore we note that the 
occurrence of a third barrier in constrained mean-field 
calculations of PES’s of actinides depends on the specific 
choice for the energy density functional. 
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PES have been obtained by examining the neutron and 
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